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ABSTRACT 

We compile a sample consisting of 56 radio-quiet active galactic nuclei so as to investigate statistical properties 
of hot corona of accretion disks from ASCA observations. The black-hole masses in the sample are estimated via 
several popular methods and the bolometric luminosities from the multi-wavelength continuum. This allows us 
to estimate the Eddington ratio (£ = LBoi/^Edd) so that the undergoing physical processes can be tested via hard 
X-ray data. We find a strong correlation between Tx = ia-iokev/^Boi and £ as Tx oc £~° M with a multivariate 
regression. This indicates that the release of gravitational energy in the hot corona is controlled by the Eddington 
ratio. On the other hand, the correlation between the hard X-ray spectral index (T) and £ depends critically on 
the types of objects: T is nearly constant (T oc £°) in broad-line Seyfert l's (BLSls), whereas T oc log£ 018 in 
narrow-line Seyfert l's (NLSls), although not very significant. We can set constraints on the forms of magnetic 
stress tensor on the condition that Tx is proportional to the fraction / of gravitational energy dissipated in the 
hot corona and that / is proportional to magnetic energy density in the disk. We find that the shear stress tensor 
f r oc P gas is favored by the correlation in the present sample, where P gas is the gas pressure. 

Subject headings: accretion, accretion disks - magnetic fields - galaxies: active - galaxies: Seyfert 



1. INTRODUCTION 

Accretion onto a supermassive black hole is generally re- 
garded to be powering active galactic nuclei (AGNs). For a 
steady Keplerian accretion disk around a black hole with mass 
Mbh, the released energy via viscosity dissipation is given by 

->25.^„,-l „-3 <r/„\ -2\ ^ 



cm 2 ) 



£ vis = 8.6 x lO^rnni^r 3 J(r), (erg s" 

where m = Mc 2 /Lem, the Eddington luminosity Lgdd = 
4-ncGM BH /K es , n es = 0.34, r = Rc 2 /2GM B h, J{r) = 1 -(3/r) 1 / 2 
and »jbh = Mbh/Mq (Shakura & Sunyaev 1973). It is impor- 
tant to note that £>vis is independent of the specific mechanism 
of viscosity stress. In the regime of the standard accretion disk 
around a massive black hole, the disks can not produce X-ray 
emission. Models of the accretion disks with hot corona have 
been suggested in order to explain the observed X-ray emis- 
sion (Liang & Price 1977). In a popular model, a fraction / 
of 2vis is assumed to be released in the hot corona and the left 
( 1 -/) in cold disk, but the transportation of angular momentum 
is mainly taking place in cold disk (Haardt & Maraschi 1991; 
Svensson & Zdziarski 1994, Kawaguchi et al. 2001). However, 
the factor / remains open since we poorly understand how the 
process happens in the disk. 

The important roles of the magnetic field turbulence have 
been realized in transportation of angular momentum (Balbus 
& Hawley 1991, Tout & Pringle 1992, Turner et al. 2002, 2003) 
and formation of the hot corona (Galeev et al. 1979, Stella & 
Rosner 1984, Merloni & Fabian 2002, Liu et al. 2002, Turner 
2004, Kuncic & Bicknell 2004). The Maxwell stress is capa- 
ble of transporting the angular momentum in the disk (Balbus 
2003). Numerical simulations show the actual physical pro- 
cesses inside the disks are very complicated (Balbus & Hawley 
1991, Stone et al. 1996; Blaes 2002; Blaes & Socrates 2003; 
Machida & Matsumoto 2003; Igumenshchev et al. 2003; Sano 
et al. 2003, Kato et al. 2004). We noted that these simulations 
are based on the assumption of radiative inefficiency, and can 
not be directly applied to explain observations. The seminal 
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paper by Galeev et al. (1979) showed that the gravitational en- 
ergy will be mainly released in the hot corona since the strong 
buoyancy and magnetic field reconnection inevitably lead to the 
formation of the hot corona above cold disk (Stella & Rosner 
1984). The formation of the hot corona is most likely naturally 
related with the transportation of angular momentum, the factor 
/ can be thus obtained if the magnetic stress and energy trans- 
portation are assumed (Merloni & Fabian 2002). This may lend 
us a possible opportunity to test the working magnetic stress 
from hard X-ray observations. There is growing evidence in the 
Galactic black hole candidates for that the hot corona becomes 
weak when m increases (see a review of McClintock & Remil- 
lard 2004). However it is elusive how the hot corona changes 
with m in AGN disks. Two problems remain open: 1) what 
is the proper magnetic field stress? 2) how to connect the two 
processes of angular momentum transportation and the energy 
release in the hot corona? 

In this Letter, we look into the hard X-ray data of radio-quiet 
AGNs to see m-dependence of hard X-rays and consider the 
roles of magnetic turbulence in the disks. We found a strong 
correlation between the factor / and the Eddington ratio, im- 
plying that the energy release is driven by the Eddington ratio 
and the importances of magnetic turbulence. 

2. THE SAMPLE AND CORRELATIONS 

Table 1 gives the data of the present sample, which is com- 
posed of 56 objects consisting of 29 narrow line Seyfert 1 galax- 
ies (NLSls) and 27 broad Seyfert 1 galaxies (BLSls). Col (1) 
is the name of the objects, Col (2) the 2-10keV luminosity from 
ASCA observations in units of erg s" 1 , Col (3) the photon index 
between 2-10 keV, Col (4) the bolometric luminosity in units 
of erg s , Col (5) the black hole mass in the units of solar 
mass. Col (6) gives some notes on the objects, indicating the 
references of X-ray luminosity and photon index as well as the 
references or methods of the estimations of the black hole mass, 
bolometric luminosity. The black hole masses are estimated by 
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Table 1 The Sample 



Name 


logLx 


r 


logL Bo i 


log JflBH 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Mrk 335 


43.42 


1.94 


44.69 


6.69 


T, 1 


IZw 1 


43.73 


2.25 


45.47 


7.26 


T, 2 


Ton SI 80 


44.00 


2.46 


45.70 


6.91 


T, 3 


F9 


44.26 


1.91 


45.23 


7.91 


T, 1 


RX J0148-27 


43.81 


2.06 


45.68 


6.97 


T, 3 


Nab 0205+024 


44.52 


2.27 


45.45 


7.86 


T, 1 


Mrk 1040 


42.83 


1.56 


44.53 


7.64 


T, 1 


LB 1727 


44.59 


1.56 


46.37 


8.47 


T, Rev 


3C 120 


44.34 


1.89 


45.34 


7.42 


T, 1 


Akn 120 


44.07 


1.91 


44.91 


8.27 


T, 1 


MCG+8-11-11 


43.54 


1.56 


44.16 


7.24 


T, Rev 


1H0707-495 


42.95 


2.27 


44.43 


6.30 


T, Rev 


NGC 3227 


42.01 


1.61 


43.86 


6.76 


T, 4 


RE 1034+396 


43.05 


2.35 


44.52 


6.45 


T, 3 


NGC 3516 


43.43 


1.83 


44.29 


7.36 


T, 1 


NGC 3783 


43.25 


1.70 


44.41 


6.94 


T, 1 


Mrk 42 


42.54 


2.14 


43.91 


5.61 


T, Rev 


NGC 4051 


41.56 


1.92 


43.56 


6.13 


T, 1 


NGC 4151 


43.01 


1.57 


43.73 


7.13 


T, 1 


Mrk 766 


43.08 


2.16 


44.23 


6.05 


T, 3 


NGC 4593 


43.06 


1.78 


44.09 


6.91 


T, 1 


IRAS 13224-3809 


43.11 


1.97 


45.74 


6.75 


T, Rev 


MCG-6-30-15 


43.11 


2.02 


43.59 


6.21 


T, Rev 


EXO 055620-3 


43.97 


1.70 


44.37 


7.21 


T, Rev 


IC 43 29 A 


43.06 


1.71 


44.78 


6.77 


T, 1 


Mrk 279 


43.99 


2.04 


44.55 


7.83 


T, Rev 


PG 1404+226 


43.49 


2.07 


45.29 


6.76 


T, Rev 


NGC 5548 


43.76 


1.79 


44.83 


8.03 


T, 1 


Mrk 478 


43.83 


2.06 


45.42 


7.27 


T, Rev 


Mrk 841 


43.83 


1.77 


45.84 


8.10 


T, 1 


Mrk 290 


43.53 


1.68 


44.35 


7.04 


T, Rev 


3C 390.3 


44.41 


1.63 


44.88 


8.55 


T, 1 


Mrk 509 


44.38 


1.82 


45.03 


7.86 


T, 1 


Akn 564 


43.38* 


2.70 


44.47 


6.04 


T, Rev 


MCG-2-58-22 


44.43 


1.73 


45.30 


8.48 


T, Rev 


E 0015+162 


45.10 


1.99 


45.31 


8.31 


R, 5 


PHL 909 


44.24 


1.11 


45.47 


8.90 


R, 1 


HE 1029-1401 


44.44 


1.81 


46.03 


9.08 


R, 1 


PG 1114+445 


44.01 


1.71 


45.65 


8.41 


R, 2 


PG 1116+215 


44.47 


2.09 


46.22 


8.50 


R, 2 


PG 1211+143 


43.70 


2.06 


45.89 


7.88 


R, 2 


PG 1216+069 


44.70 


1.57 


46.45 


9.17 


R, 2 


PG 1404+226 


42.95 


1.77 


45.27 


6.74 


R, Rev, 6 


PG 1416-129 


44.48 


1.78 


45.90 


8.50 


R, 2 


PHL 1092 


44.10 


1.99 


46.07 


7.61 


V, Rev 


RX J0439-45 


44.04 


2.25 


45.66 


7.12 


V, Rev 


PKS 0558-504 


44.67 


1.81 


45.80 


7.56 


V, Rev 


PG 1244+026 


42.10 


2.35 


44.62 


6.29 


V, 2 


PG 1543+489 


44.16 


2.46 


46.11 


7.79 


V, Rev 


IRAS 1702+454 


43.74 


2.20 


44.74 


6.54 


V,Rev 


Mrk 507 


42.62 


1.61 


44.54 


6.34 


V,Rev 


IRAS 20181-224 


43.87 


2.33 


45.43 


6.80 


V, Rev,7 


Mrk 142 


43.18 


2.12 


44.55 


6.65 


L, 3 


IRAS 13349+2438 


44.12 


2.31 


45.35 


7.74 


L, 3 


Kaz 163 


43.19 


1.92 


45.02 


7.11 


L, Rev 


NGC 7469 


43.26 


1.91 


45.28 


6.84 


N, 1 



NOTES/REFERENCES:- T: Turner et al. (1999); R: Reeves & Turner (2000); 
V: Vaughan et al. (1999); L: Leighly (1999); N: Nandra et al. (2000); 1: Woo 
& Urry (2002); 2: Vestergaard (2002); 3: Grape et al. (2004); 4: McLure 
& Dunlop (2001); 5: Orndahl et al. (2003); 6: Boroson & Green (1992); 7: 
Elizalde & Steiner (1994); Rev: reverberation relation from the magnitude Mb 
from Veron-Cetty & Veron Catalog; * data is taken from Turner et al (2001). 

1) the empirical reverberation relation with calibrations, Rblr = 
30.2 (L 5 i 00 /10 44 ) lt-days, where L 5 ioo is the luminosity at 
5 100A( Vestergaard 2002); 2) reverberation mapping technique 
(Kaspi et al. 2000); 3) relation between BH mass and its host 
galaxy magnitude logM Bn /M Q = -0.5M R -2.96 (McLure & 
Dunlop 2002); 4) relation between the dispersion velocity and 
BH mass M BH = 1.35 x 10V/200) 4 ' 02 (Tremaine et al. 2002). 
The black hole masses in Table 1 are taken from Woo & Urry 
(2002) and Vestergaard (2002), otherwise we use the empirical 
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FIG . 1 . — The plots of hard X-ray spectrum index and the ratio of ASCA lumi- 
nosity to bolometric luminosity versus the Eddington ratio. The open and filled 
circles represent the narrow line Seyfert 1 galaxies (FWHM H/3 < 2000km/s) 
and the broad Seyfert 1 galaxies, respectively. The object IRAS 13324-3809 
is a NLS1 with FWHM(H/3) = 620km/s and 2-10keV luminosity 10 43 erg/s 
much fainter than its Mb = -24.2 (Veron-Cetty & Veron's catalog), with 
m BH = 10 6 ' 95 , the ratio Lj-iokeV/^Edd ~ 10~ 3 and L Bo i/L E dd ~ 10° 88 . 

reverberation method to estimate Mbh- £5100 is obtained from 
the absolute M B given in Veron-Cetty & Veron's quasar catalog 
(1 1th edition) through an extrapolation of a power law spectrum 
in optical band as F v oc 5 . The bolometric luminosities LboI 
are taken from Woo & Urry (2002) and Grupe et al. (2003), oth- 
erwise they are estimated via L Bo i = 9L5100 (Kaspi et al. 2000). 
We use the Hubble constant Ho = 75km s" 1 Mpc" 1 and deceler- 
ation factor qo = 0.5 in this paper. 

Fig 1 shows the plot of photon index, T, and fraction of 
X-ray emission, T-g. = i^-iokev/^Bob versus the Eddington ra- 
tio, £ = LBoi/^Edd- Except for PHL 909, the X-ray photon in- 
dices almost keep a constant for BLSls with a mean value of 
(r) = 1.78±0.17. There is a trend for NLS Is: the higher is ac- 
cretion rate, the steeper is the X-ray spectrum as T oc log£ 018 . 
We find the Pearson's correlation coefficient of T — £ is 0.37 
for NLS Is. The difference of £ — T relations in BLSls and 
NLS Is show an interesting evidence for the bi-modal accre- 
tion in Seyfert 1 galaxies. However, it is interesting to note that 
there is a strong correlation of T and the Eddington ratio £ for 
the entire sample, 

T = (2.05 ± 0.04) + (0.26 ± 0.05) log£, (2) 

with the Pearson's coefficient r = 0.61 and probability p = 
5.4 x 10~ 7 . We use the least square method of multivariate 
regression and find a very strong correlation logLx = (0.36 ± 
0.09) logLfioi + (0.47 ± 0.09) logL Edd + (6.13 ±3 .57) with coef- 
ficient r = 0.83 and F = 58.7 of the F-test (much larger than 



Hot Disk Corona and Magnetic Turbulence 3 



Table 2. Magnetic Turbulence as Viscosity and Properties of Hot Corona 



viscosity stress 


Re fere nc e s 


factor f 


P i ^ P 

*rad ^ -'gas 


P 

fga 


^ P i 


Comments 


1) tr(j> = -OtPvn 


SS73 


ill 


/ = const. 


/ = 


const. 


X 


2) t r( p = -afgas 


SR84 


oc a l l 2 a s /c s 


/ oc m" 7 


/" 


•» const. 




3) t r4> = -aP Ia d 


WL91 


cxaV2[l_(a s / Cs fl l ' 2 


fx a 1 / 2 


/" 


■» 


X 


4) t r<p = -aV-Pradfiot 


S90 


cxaVzfl-fe/c^]'/" 


/oca'/2 


/" 


■» 


X 


5) t r( p = —0<-\J /gasAot 


TL84, B85, MF02 


ex a'^iyoj/c, 


/ oc rir 1 ! 1 


/" 


-» const. 


X 


6)f r0 =-a V / Pga S / , rad 


LN89 


X a 1 ^a./c\l-(fl,/e,f\ 1/4 


f oc a l l 2 m~tl 2 


/- 


^0 


X 



REFERENCES:- B85: Burm (1985); MF02: Merloni & Fabian (2002). LN89: Laor & Netzer (1989); S90: Szuszkiewicz (1990); 
SS73: Shakura & Sunyaev (1973); SR84: Stella & Rosner (1984); TL84: Taam & Lin (1984); WL92: Wandel & Liang (1992). 
Here 7 is a positive constant, and a s = (Pgas/p) 1 ^ 2 the gas sound speed. 



F = 3.2 at 99.9% level of confidence). This correlation can be 
rewritten as 



logJSc = -0.641og£-0.171ogL Edd + 6.13, 



(3) 



from which an interesting conclusion follows: the hot corona is 
mainly controlled by the Eddington ratio. Fig lb. shows this 
correlation. 

The correlation shows that the corona becomes dramatically 
weak with increases of the Eddington ratio. While the Edding- 
ton ratio is low, a larger fraction of the energy will be released in 
the hot corona. A naive expectation is that the magnetic stress 
responsible for the transportation of angular momentum in- 
creases with the accretion rates. Otherwise the accretion would 
be halted because of the inefficient transportation. However, 
this does not agree with the observed correlation, which indi- 
cates that the dissipation in the cold disk (not in the hot corona) 
becomes more efficient with increases of the accretion rates. 
Successful models of the hot corona should address this issue 
explicitly. 

The correlation of £ - T shows that the hard X-ray spectrum 
becomes harder with increases of £, confirming the results of 
Lu & Yu (1999) and Gierlinski & Done (2004). This is usu- 
ally explained by that the seed photons efficiently cool the hot 
corona with increases of m. However, this explanation could 
be replaced by that the hot corona really becomes weak since 
^"x oc £-° M . Indeed this agrees with that the broad line Seyfert 
1 galaxies have strong reprocessing of hard X-ray, such as in 
NGC 5548 (Clavel et al. 1992; Chiang et al. 2000), evidenced 
by that optical and UV variations follow the hard X-ray. In 
NSLls the reprocessing is rather weak, such as Akn 564 (Shem- 
mer et al. 2001), NGC 4051 (Done et al. 1990) and IRAS 
13224-3809 (Young etal. 1999). 



3. MAXWELL STRESS AND HOT CORONA 

Assuming that the energy transportation from the cold disk 
to the hot corona is by buoyancy driving the magnetic tubes 
across extent of the disk, Merloni & Fabian (2002) proposed 
the fraction of the energy transported by magnetic buoyancy to 
be / = P mll gV p /Q+, where v p is the transporting velocity, P mag 



the magnetic pressure and the dissipated energy Q+ 
via viscosity stress t r< p. We have 



2b v A 
3ko c s 3ko 



2 3 / 2 b / P, 



1/2 



(4) 



where f,.^ = -koP mllg , = B/(47rp) 1//2 Alfven velocity and the 
constant b = v p /va- Here P tot = P ra d + P ga s and c 2 s = P tol / p. 
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FIG. 2. — The plot of factor / and the Eddington ratio in theoretical models. 
The lines labeled by MF02 and SR84 are results of Merloni & Fabian (2002) 
and Stella & Rosner (1984), respectively. We take 2 3 / 2 b/3k = 1 and b = 1. 
The solid, dashed, and dot-dash-dot-dash lines are results for a = 0.05,0.1,0.4 
for bibh = 10 8 . The two thick lines are for mgH = 10 6 and a = 0.05. £ = rjm is 
used, where the accretion efficiency r\ — 0.1. 

According to eq.(4), Table 2 gives the properties of the factor 
/ for six distinct types of the magnetic stresses based on differ- 
ent assumptions. Comments are given in Table 2 based on the 
correlation (eq. 3). Cases 1), 3), 4) and 6) predict properties in- 
consistent with the correlations in Figure lb, thus they could be 
ruled out. Only Cases 2) and 5) are possible since they provide 
characters of the hot corona consistent with the correlation. 

Specifying the magnetic stress, we can solve the structure of 
disk and the factor / at each radius formulated in Svensson & 
Zdziarski (1994) as described by Merloni & Fabian (2002) and 
Merloni (2003). We define the radial-averaged fraction (/), 

^ f(r)Q + (r)rdr 



(f) 



f™Q + (r)rdr 



(5) 



so that we can compare with the observation, if the magnetic 
stress is assumed. 

In Case 2) the magnetic pressure t r< p = -aP gas is originally 
suggested by Stella & Rosner (1984), who argued the buoyancy 
controls the growth of the magnetic field. The magnetic stress 
is simply scaled by the gas pressure. Such a stress can stabilize 
the thermal instability of the cold disk (Svensson & Zdziarski 
1994). On the other hand, Ichimaru (1975) argued that the 
growth of the magnetic field may be controlled by the recon- 
nection and gave a magnetic pressure tensor t r< f, = -a^/P gas P tot 
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(Case 5). This stress has been extensively studied by many au- 
thors (Burm 1985, Szuszkiewicz 1990). Recently Merloni & 
Fabian (2002) and Merloni (2003) derived this simple formula- 
tion based on analytical work (Blaes & Socrates 2001) and 3-D 
MHD numerical simulation (Turner et al. 2002) of magneto- 
rotational-instability in radiative accretion disks. 

Figure 2. shows the results of (/) for the magnetic stresses 
t r< p = -afgas and t r( p = -a^/P gas P tot . We find that (/) is not sen- 
sitive to the black hole mass. More interesting thing in Figure 2 
is that the slope of £ — (f) (for £ > 0.05) is insensitive to both 
the black hole mass and the viscosity a, but only sensitive to the 
magnetic stress. These properties of the theoretical relation of 
£— (f) allow us to test the working magnetic stress in disk. We 
have (/) oc £~ - 44 for the Merloni & Fabian's model whereas 
(/} oc £~ QJ1 for the Stella & Rosner's. MF02 predicts a flatter 
slope than the observation. SR84 model has a steeper slope, but 
it matches the observational correlation (/) oc £-°- 64 ± 009 within 
the uncertainties in the present sample. 

A more fiducial test on the magnetic stress requires the accu- 
rate estimations of the black hole masses, bolometric luminos- 
ity and the total hard X-ray luminosity from observations. The 



near future observations of INTEGRAL can provide entire hard 
X-ray spectra of a large sample so as to test the disk-corona 
models. The future 3-D MHD numerical simulations including 
the hot corona should be tested by the plot of £-L x /Z-Boi- 

4. CONCLUSIONS 

With the estimations of the black hole masses in radio-quiet 
AGNs, we find the hard X-ray spectrum index and Z^-iokev/^Boi 
strongly correlate with the Eddington ratio. These correla- 
tions directly indicate that the fraction / is controlled by the 
Eddington ratio. Such a correlation can be explained by the 
magnetic turbulence that plays two roles in transporting angu- 
lar momentum and forming the hot corona. The correlation 
Tx oc £-°- 64 ± - 09 favors the magnetic stress t r< f, oc P gas in the 
present sample. 

This research is supported by Grant for Distinguished Young 
Scientists from NSFC, 973 project and a Grant-in-Aid of the 
21st Century COE "Center for Diversity and Universality in 
Physics". 
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